The Energy- Weighted and Non Energy- Weighted Gamow- Teller Sum Rules in 
Relativistic Random Phase Approximation 
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The non energy-weighted Gamow- Teller(GT) sum rule is satisfied in relativistic models, when 
all nuclear density-dependent terms, including Pauli blocking terms from nucleon-antinucleon ex- 
citations, are taken into account in the RPA correlation function. The no-sea approximation is 
equivalent to this approximation for the giant GT resonance state and satisfies the sum rule, but 
each of the total 0- and j3+ strengths is different in the two approximations. It is also shown that 
the energy-weighted sum of the GT strengths for the f3- and f3+ transitions in RPA is equal to the 
expectation value of the double commutator of the nuclear Hamiltonian with the GT operator, when 
the expectation value is calculated with the ground state in the mean field approximation. Since 
the present RPA neglects renormalization of the divergence, however, the energy-weighted strengths 
outside of the giant GT resonance region become negative. These facts are shown by calculating in 
an analytic way the GT strengths of nuclear matter. 

PACS numbers: 21.60.-n, 21.60.Jz, 21.65.+f 



I. INTRODUCTION 

Recently, the present authors have investigated Ikeda- 
Fujii-Fujita(IFF) sum ruleQ of Gamow-Teller(GT) tran- 
sitions in relativistic models^, It has been shown 
that the sum of the GT strengths in the nucleon sec- 
tor is quenched by about 6% in finite nuclei^] and 12% 
in nuclear matter 0- This fact is owing to the small 
component of nucleon wave functions, and the quenched 
amount is taken by nucleon-antinucleon states. Although 
the nucleon-antinucleon states are far from the giant GT 
resonance region, if there is a coupling with particle- hole 
states, they may contribute to nuclear excitations in the 
giant resonance region as virtual states. In the previous 
paper 0, 0, the effects of the coupling were estimated 
with use of the random phase approximation(RPA) in rel- 
ativistic nuclear models. RPA correlation functions were 
calculated by neglecting the nuclear density-independent 
terms in order to avoid the divergence, but by keeping all 
the density-dependent terms including the Pauli blocking 
ones from nucleon-antinucleon excitations. This method 
was frequently used in the study of no charge- exchange 
excitations |^,|5|,|6|. From now on let us call this approxi- 
mation no free term approximation(NFA). In this calcu- 
lation, we found that effects of the coupling are negligible 
on the excitation energy and strength of the GT state in 
the nucleon sector. 

The purpose of the present paper is threefold. First, 
we will study whether or not NFA satisfies IFF sum rule 
which is the non energy-weighted sum rule for the differ- 
ence between the /?_ and /3+ transition strengths. When 
we discuss excitation strengths, the used approximation 
should satisfy at least the model-independent sum rule. 
It seems not to be trivial for NFA to satisfy IFF sum 
rule, since for the /?_ transitions, the backward ampli- 
tudes in the particle-hole sector do not contribute to the 



strength as in Tamm-Dankoff approximation, while in the 
nucleon-antinucleon excitations, they are comparable to 
the forward ones in the present model. We will show, 
however, that the sum rule is satisfied in NFA. The proof 
will be done in an analytic way for nuclear matter. 

Second, we will discuss the sum rule in the no-sea ap- 
proximation(NSA), which is also used frequently in no 
charge-exchange excitations 0, HI- It will be shown that 
NSA also satisfies IFF sum rule, but that each total 
strength of the /3_ and /3 + transitions is different from 
those of NFA. As far as the giant GT resonance state, 
however, its excitation energy and strength are predicted 
in the same way as in NFA. 

Third, we will investigate the energy-weighted sum 
of the GT strengths in RPA, which have not been dis- 
cussed so far, as far as the authors know. For no charge- 
exchange excitations, there is a famous theorem in the 
non-relativistic framework that the energy- weighted sum 
of the excitation strengths in RPA is equal to the ex- 
pectation value of the double commutator of the nu- 
clear Hamiltonian with the relevant operator 0. Here 
the expectation value is calculated with the Hartree-Fock 
ground state. It will be shown in relativistic models that 
the same theorem holds for charge-exchange excitations 
with respect to the sum of the GT strengths for /3_ and 
j3 + transitions. In both NFA and NSA, however, the 
energy-weighted strength outside of the giant GT reso- 
nance region is negative. This is due to the fact that 
divergent terms are simply neglected in the two approxi- 
mations. We definitely need the renormalization to solve 
this problem in future work. 

We note here that such detailed discussions as men- 
tioned in the above are required not only from a theoret- 
ical point of view, but also from the recent experiment. It 
has been shown experimentally that the sum rule value is 
quenched by about 10% from the non-relativistic analysis 
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of (p,n) reaction^(j- So far, all the 10% quenching was 
assumed to be due to the coupling of the particle-hole 
states with A-hole states 0, u^- If the 6% quenching 
is due to the relativistic effects, it may conclude that 
only 4% stems from the contribution of the A degrees 
of freedom. Such a weak coupling of the particle-hole 
states with the A-hole states may yield the extremely 
small value of Landau-Migdal parameter <?^ A which dom- 
inates the critical density of the pion condensation^ 
and other spin-dependent response functions 0]. Thus 
we need to discuss GT strengths carefully, before study- 
ing spin-dependent structure of nuclei in detail. 

In the next section, we will briefly review NFA and 
NSA in no charge-exchange excitations. In section IIIII 
the mean field correlation function of the GT excitations 
will be derived. In section Hvl we will discuss the GT sum 
rules in RPA within the framework of NFA and NSA. 
The final section will be devoted to discussions and con- 
clusions of the present paper. 



II. NFA AND NSA 

The fundamental problem in relativistic nuclear mod- 
els is how to renormalize the divergence due to the anti- 
nucleon degrees of freedom. Only few attempts of the 
renormalization have been reported so far 0,^3- Never- 
theless, it has been shown for the past 30 years that rela- 
tivistic models explain phenomenologically very well nu- 
clear structure and reactions without renormalization |lfj 
Il7| . In those studies, there are some cases where it is 
necessary to take into account a part of the antinucleon 
degrees of freedom. In RPA based on the mean field ap- 
proximation, the continuity equation of the baryon cur- 
rent is violated, if the configuration space is limited to 
the nucleon sector only. At present, there are two ways 
to avoid the violation. The one(NFA) is to neglect the 
density-independent terms in the RPA correlation func- 
tion which are divergent, but to keep all the density- 
dependent terms including the Pauli blocking terms from 
nucleon-antinucleon excitations |J, |5|, |6j- The other is 
called the no-sea approximation(NSA), where the Dirac 
sea is assumed to be empty, and the antinucleon states 
are treated as particle states with negative energies in 
the configuration space of RPA0. NSA is described by 
changing the sign of the imaginary part of the Green 
function of the antinucleon. It is shown in this way that 
the divergence of the RPA correlation function disap- 
pears and that the continuity equation is not violated^- 

In this section, we briefly review the above two ap- 



proximations in nuclear matter, since their structure in 
charge-exchange excitations is a little different from that 
in no charge-exchange ones discussed so far. 

The Green function of the Dirac particle in nuclear 
matter is written asP.H^|. 

G(p) = (1 - P ) G p (p) + 9 P G h ( P ) + G$(p), (1) 

where the particle, hole and antinucleon parts are given, 
respectively, as 



G P (P) 
G h (p) 



^+(P) 



Gn(p) = - 



Po - E p + ie 
Po - E p - ie 



Po + E p - ie 



(2) 
(3) 
(4) 



The projection operators in the above equation are de- 
fined as usual, 



^+(P) = ^2u a (p)u a (p) 



f + M* 



2E n 



A_(p) = -^u a (-p)u a (-p) 



p + M* 
2E 



with 



p» = (E p ,p), r = {~E P ,p), E p = v/M*2+p2. 

The step function is expressed by using the abbreviation: 
6p = 9(k F — |p|), fcp being the Fermi momentum. The 
Green function Eq.Q) can be also expressed in terms of 
the density-dependent and -independent parts as, 

G(p) = G D ( P ) + G F (P), (5) 
G D ( P ) = 6 p (G h (p)-G p (p)) , (6) 
G F (P) = G p (p) + Gu(p). (7) 

The mean field correlation function is described as 
1 



n 



2ni 



d^pTv{ la G{p + q ) lb G{p)), (8) 



where 7 and jt> stand for the 4x4 matrix of the external 
field. If Eq.© is inserted into Eq.^SJ), the mean field 
correlation function is divergent, because of the density- 
independent term which contains GfG f . In NFA, it is 
simply neglected without renormalization, to have 



iT D = - ^ / d 4 P Tr la (G D (p + qhbG B (p) + G F (p + qhbG u (p) + G B (p + q) lb G F {p)) . (9) 
Keeping the terms which remain after integration over po in the complex plane, the above equation is separated into 
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two parts, 

77 D = 77 P h + TTpauii. (10) 
The first term represents the particle-hole correlation function, while the second term the nucleon-antinucleon one, 

tfph = - ^- J d 4 pTr la ((1 - p+q ) P G p (p + q) lb G h { P ) + (1 - P ) fl p+q G h (p + q) lb G p {p)) , (11) 

^Pauii = ^.[d A pTr la (e p+q G p (p + q) lb G$(p)+d p G$(p + q)j b G p (p)). (12) 

Thus, ilpauii is composed of the Pauli blocking terms. In NFA, however, these terms are necessary for keeping the 
continuity equation of the baryon current, and for satisfying the GT sum rule, as shown later. 

In NSA, Green function of the antinuclcon in Eq.(0J) is artificially modified by changing a sign of the imaginary 
part, 

Gbo(p) = - (13) 
P o + E p + is 

Then density- independent terms have no contribution to the integration over P q in Eq.JSJ, and we can avoid the 
divergence problem without violating the continuity equation^ ■ The terms which should be kept in NSA are written 
as 

n = 77 ph + i7 no „ sca , (14) 
where the particle-hole correlation function is the same as in NFA, but the antinucleon-dependent part is given by 

-ffno-sca = - tt~ d 4 P Tr~f a (6> p+q G h (p + q)~f b G no (p) + 6 p G no ( P + q)^ b G h (p)) . (15) 
2m J 

After integration over P q, the mean field correlation functions in NFA and NSA are finally written, respectively, as 

n Ph = /^Tr (7 ^ + (p + q ) 7 ^ + ( p)) ( - d-jWflp . - d-gp)gp +q ) j (16) 

J \E p +q - E p - q a - is £/ p+q - E p - q + is) 



IdpnuU = J d 3 P Tr ^6< p+q 



'p+q — -c/p — qa — «t ^p+q — & p 
7 a yl + (p + q)7 b yl_(p) 7a^-(p + q)7b^+(p) ' 



-E p+q + E p - q - ie p E p+ q + E p + q Q - ie J ' 



(17) 



77 - / J 3 nTr I ft 7a^+(P + q)7&^-(p) , „ la A- (p + q)7b^+ (p) \ , ^ 

iino-sea - / <t P it #p+q ■ + #p — — . (18) 

^p+q + -£%> — 9o + 2£ ^p+q + E p + qo + «e , 



r 



Thus, the real parts are the same in NFA and NSA, while 
the imaginary parts of -Z7p au i; and i7 no -sca have an oppo- 
site sign to each other. This fact implies that both ap- 
proximations yield the same excitation energies for the 
discrete states in RPA, but there is a possibility that 
their response functions and sum values of the excitation 
strengths are different from each other. An example will 
be shown later in the case of the GT strengths. 

We note that the relationship between the Landau- 
Migdal parameters and the correlation functions is the 
same in NFA and NSA, since it depends on the only real 
part of the correlation functions, as shown in ref . [a. Il8|| . 
All physical quantities expressed in terms of the Landau- 
Migdal parameters, therefore, must be the same in the 
two approximations. The reason why NSA as well as 
NFA does not violate the continuity equation is because 



it is also independent of the imaginary part. This fact is 
verified by showing explicitly 

# P h + ilpauli = Uph + ^no-sea = (19) 

for the operator 7b = 4- Each of them for 7/, = ^ is 
written, respectively, asQ 

n * h = J ^ 4£ p Ep+q ^ (7a7 ° ^ " #/)) ' 

-^Pauli — -^no-sea 

= J d3 P AE J ~ Tr (7a 7o (0 - , 

using the notations: 

j/" = (£ p+q , p + q), = (-£ P +q, p + q). 



III. THE MEAN FIELD CORRELATION 
FUNCTIONS FOR GT EXCITATIONS 

When we discuss GT excitations in N ^ Z nuclei, we 
need the following replacement in Eas. (|16|) . (|17|) and (I18II . 

la -> 7a = Ibly T+ , Jb ~> 7b = 757y T - ) 

61 _> 1 Z Tz g(p) i 1 ± Tz g(°) 
p p ^ 2 p 2 p ' 

where we have defined the isospin operators for the con- 
venience as 

t± = (t x ± tYj,) /\/2 , t = r 2 , 

and the step function with respect to the proton(fc p ) and 
neutron Fermi momentum(fc„) : 

9® = 6(ki-\p\), i=p,n. 

Then, for the /3_ excitation at q — in N > Z nuclei, 
Ea. ((TB)l becomes to be 



77 ph = -4 jd 3 p 



M* 2 +pgflW _g(p) 

£2 g +i£ ' 



(20) 



and Eas. ((T7|l and (|T%|) are described, respectively, as 



/p 2 — r> 2 
d 3 P E2 y 

( of 



qo - 2E P + ie qo + 2E P - ie 




11 - - I I d 3 v p2 ^ 

JJ no-sca — ^± I u, jj 



El 



9 (n) 



(P) \ 

p 1 (22) 



V qo + 275 p + ie g - 275 p - is 

From the above equations, we see, at this stage, first that 
the present model has only the forward amplitudes for 
particle-hole pairs. Second, in NFA, the forward ampli- 
tudes of the nucleon-antinucleon pairs have an opposite 
sign to the one of the particle-hole pairs. This means 
that the GT strength from the Pauli blocking terms is 
negative. Third, in NSA, the excitation energy of the 
antinucleon-hole pairs is negative, as mentioned before. 
Finally we can see that the density dependence of NSA, 
expressed by the step functions, is different from the one 
of NFA. As a result, their GT strengths are different from 
each other, although the difference between the strengths 
of (3- and /3+ transitions is the same in the two approxi- 
mations. 

The correlation functions for the /3+ transitions are 
obtained by changing the sign of go m Eas. i|2(J[l to l]22p. 
In this case, there are no forward amplitudes for particle- 
hole pairs, but backward ones. 



By performing the integration, Eg. (|20fl is written as 

n P h = —7— , a p h = (Q(K) ~ Q{k p )) , 

qo + i£ 3 



where Q(ki) is defined by 



An 



d 3 p 



M* 2 + p 2 y 

E 2 
p 



-i ki 



+ 2Mf - 2M * 3 tan — -■ (23) 
3 M* 



Eqs.J2U an( i l|22|l are expressed by separating into the 
real and imaginary parts: 



77 F 



77, 



(R) 



hi, 



(i) 



iTT ( 

no-sca 1 no-scaJ 



Pauli 
(I) 



'Pauli — -'-'Pauli 

77 - 77 (R) 

"no-sca — no-s 

which are related with each other as 

7t(R) = 7t(R) n m = _TT(I) 

Pauli no-sca ' "Pauli "no-sea- 

The explicit forms of 77p R uli and i7p2uli are described as 



77 



(R) 

Pauli 



167T 



P N( k p> - ?o) + P^(kn, qo)) , 



-^Pauli = 7 Pau ii ( fc p , g ) + 7 Pau ii ( fc n , - g ) , 

where we have defined the two functions: 



Pjr(k F ,q ) = — 



' d 3 p P 2 - p 



An J El 2E p + q ' 
7p au ii(fcF,go) = -An J d 3 p E2 y 5{q - 2E p ) 

The second function is simply expressed as 
7p au ii(fcF, go) 

r*2\ 3 / 2 



4tt 2 (gg - AM* 
3 go 



0. 



2M* < g < 2E F , 
otherwise, 



(24) 



with Ep = \J M* 2 + fcp, while the first one is written as 
P^{kF,qo) 



' F 

k F E F g fc F 6A7* 2 - - 2 



gg fc F + Ef 
— log- 



2M* 3 _-, fc F 
tan — — 

go M* 



4g 



A7* 



7b, 



where 7b is given by 



i-M* 



q 2 - AM* 



7r = / dt ■ „ 

J EF -k F t 2 + qot + M* 2 



- 4MV1 - x 2 tan" 1 y , x 2 < 1 
2M*\/x 2 - Hog ) ,-- 1 



1-1/ 



with 



V 



Ep — M* y/|l-a; 2 
kp 1 + x 

For g <C 2M*, we have 



go 
2M* 



:-! /" fcF d 3 p p 2 - p 2 y 



4* Jo El 2E p 
El ( -vp - vp - - (l - u|) log y 



1 +_«F 

1>F 



.2 4 



1 
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where vp denotes the Fermi velocity kp/Ep. Thus, in 
contrast to the imaginary part, the real part of 7Tp au ii and 
^7no-sca may yield a small contribution to the excitation 
energy and strength of the giant GT resonance state in 
the region g <C 2M* . 

The response function i?(go) is defined with the above 
correlation functions as[5j 



R(qo) 



l v 



7T (27r) 



Im n(q ) 



V standing for the volume of the system (3tt 2 /2kp)A. 
Then, the total GT strength in the mean field approx- 
imation is obtained by integrating it over qq. The con- 
tribution from the particle-hole states to the strength of 
the /3_ transitions is 

11^ /*°° 

s * = n^r^oj_ v dqolmIIph{qo) 



W) 



, , ,\ M* + p 



El 



> (25) 



where the variable ry is introduced, since Im77 p h(go) con- 
tains 8(qo). This can be written with use of Eq.J^BJ) as 



(-) 
P h 



A 



TJ (Q(fcn) - Q(fcp)) 



The function Q(k n ) — Q(k p ) can be expanded in terms of 
Q(fc„) - Q(fcp) « ^^(fc„ - fc p ) + • • • , 



where we have 

dQ(kp) 



dkp 



= 3ki 1 - -trf 



When utilizing as usual the relationship: 

k n ~ kp ~ ~^ ' 

the GT strength in the nucleon sector is approximately 
written as 



In the present definition of the spin-isospin operators, 
the well-known IFF sum rule is described as 



\Q+Q-\ )-< \Q-Q+\ ) =2{N-Z) (27) 



with 



i 

This sum rule is nothing but a result of the the commu- 
tation relation: [r+ciy , T-(X V ] = 2t z . If we assume that 
Q+\ ) = 0, then we simply obtain 

(\Q+Q-\)=2(N-Z). (28) 

Comparing Ea. 1)26(1 with the above equation, it is seen 
that the GT strength in the nucleon sector of the rel- 
ativistic model is quenched by a factor (1 — §«§). In 
fact, as shown below, the quenched amount is taken by 
the nucleon-antinucleon pair excitations. In the present 
model, there is no GT strength of the f3+ transitions in 
the particle-hole sector, 

1 V f°° 
S pb = - 7^7)3 ,5? J d 1o ImiI P h(-go) = 0. 

The total GT strength of the proton-antineutron(/3_) 
excitations in NFA is given by 

1 V f°° 

dq Imn Pauii (q ) 



j(-) 

'Pauli 



7T (27T)3 J 

W 
~(2^)3 



d 3 p 9^ 



K 



E l 



(29) 



which is negative, as mentioned before. The one of 
the neutron-antiproton(/3 + ) excitations is obtained in the 
same way by changing the sign of go m Eq. H21|l , 



5, 



(+) 

Pauli 



1 V f°° 

_ 7TT7 / rf ?oImiTpauii(-go) 
7T (2iry J 



W 
(27p 



d 3 p 6> p n) 



P% 



El 



(30) 



From Eqs.JHJ), JU^I and PU|). we obtain the GT sum 
value in the mean field approximation of NFA, 



^ P h ) +4;uii-4tun = 2(iV-Z), 



(31) 



which is just IFF sum rule. Thus the quenched amount 
of the GT strength in the nucleon sector is taken by the 
nucleon-antinucleon pairs, and the Pauli blocking terms 
are necessary for NFA to satisfy the sum rule. 

In NSA, the GT strength of the particle-hole pairs is 
the same as in NFA, but the strengths from the antinu- 
cleon degrees of freedom are different. For the (3- transi- 
tions, the integration of Ea. i|22|) over negative excitation 
energy gives 



o(-) 



1 V 



7T (27r) S 



dq Im77no- sea (go) 



6 



This is equal to — 5p^ ; as seen in Eq. (|3*0jl . 



c(-) 



-5' 



(+) 

Pauli ' 



In the same way, we obtain the relationship: 



<?(+) 



-5 



(-) 

Pauli' 



where the l.h.s denotes the total strength of the /3+ 
antinuclcon-hole excitations in NSA. In NSA, each 
strength is positive, but the energy-weighted sum be- 
comes negative. Although the strengths of the /3_ and 
(3 + transitions in NSA are different from those in NFA, 
their difference satisfies the sum rule, as in Ea. l|31|l . 



sl h ) +SLl a -SiH a = 2(N-Z). 



(32) 



where 



Although it is not unique how to introduce g' in the rela- 
tivistic model, we have shown that the above Lagrangian 
yields the known expression for the excitation energy of 
the GT state in the non-relativistic limit |2j, y| . We note 
that, for example, the GT state can not be described in 
relativistic models by putting the LM parameter into the 
meson propagators |3j. 

For the Lagrangian Eci. (|33[) . the RPA correlation func- 
tion TTrpa is written in terms of the mean field one U 
as. 



IV. THE GT SUM RULE IN RELATIVISTIC 
RPA 

In non-relativistic models, RPA correlations in the GT 
states are assumed to be induced by Landau-Migdal(LM) 
parameter q'flfj [2fT|. In the relativistic model, we also 
introduce g' which is provided in the nuclear Lagrangian 
as a contact term with the pseudo vector coupling: 



rr / \ n(q ) 55 
iiRPAWoJ - 7- 777 — r , X5 - TTTvi 

l + X5^7(go) (2tt) j 



Then the RPA response functions for the /?_ and (3^ 
transitions are given, respectively, by 



r kpa (lo ) = \ 7^~)3 Inl Urpa ( ±( ?o ) ■ 



(33) In NFA, II (q ) is provided by Eqs.JU and |23} as, 



n(qo 



4 d 3 p 



? 9 

p -pi 



E 2 



q + i£ 

while the one in NSA is given by Eas. (|2TH) and (|2*2*|l as 



n(qo 



qo + i£ 



4 d a p 



p 2 -pI 



E 2 
p 



q(p) 



q(p) 
7 P 



+ 



9 (n) 
7 P 



2E P -qo-ie 2E P + q Q - is I ' 



o(») 
7 p 



2i? p - qo + is 2E P + q + is 



(34) 



(35) 



A. The Non Energy- Weighted GT Sum Rule 

In NFA, the non energy-weighted GT sum value in RPA is given by 



f f-OQ />00 

Srpa = lim n ^ / d ?o -Rrpa(?o) - / dq R 
■n^+o U-,, Jr, 



(+) 

RPA 



(90) 



In order to calculate the r.h.s., we expand ZZrpa in terms of xs> 



n=0 



where In^ is defined as 



/OO />OC 
rfg Im(7T( (Z o)r +1 - / dgo Im (77 
-7} J 7] 



(-*>)) 



n+l 



(36) 



(37) 
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Since 77(±gn) has no pole in the first quadrant of the complex plane, the integration on the closed contour provides 
us withEl 



dq (77(g )f +1 - i 



dq (II(iq - V )) n+1 + I dq (77(g ))" +1 = 0, 



\n+i 



c 



dq (77(-g )) 



n+1 



dq (n(-iq - V )) n + 1 + I dq (77(-g )) 



N.n+1 



(38) 
(39) 



where C indicates the contour on go = Re 10 , (0 < 9 < tt/2). Using the above two equations together with the fact 
that Il*(—iqo — if) — 77(igo — if), 7„ can be expressed as 



7(0) 

When |g | — > oo, 77(g ) in Ea.l|34|l behaves as 

n(qo) 

Hence, the integration on the contour C gives 



Im / dq (n(q )) n+1 -(n(-q )) n+1 



2(2tt) 3 TV- 



0.0 



V 



Jc 



dq (77(±g )) 



n+l 



I 



T2(2tt) 



,N - Z 
V 



n+l r ir/2 



d6e 



0, (R^oo, n> 1 



(40) 



(41) 



(42) 



r 



From Eas. and finally we obtain 

^=0, (n>l). 

Hence, the GT sum value of RPA in NFA is equal to the 
one in the mean field approximation Ea. l|31() . 



Srpa = — 



lira /, 



(o) 



7T (27r) 3 



2(N-Z). 



Thus, RPA in NFA satisfies IFF sum rule. 

We can calculate in the same way the GT sum value 
of RPA in NSA, but by defining J^ 0) instead of ^ 0) in 



/( 



o) 



dq Im (77 (g )) 



n+l 



dq Q Im(77(-g )) 



n+l 



where 77(go) is given by Eas. (|35|) . By taking a closed 
contour with no pole in the second quadrant in this case, 
we can prove that J„ 0) = (n > 1), and therefore, NSA 
also satisfies IFF sum rule. 

In the above, we have shown that both NFA and NSA 
satisfy IFF sum rule. This means that the difference be- 
tween total GT strengths for the /3_ and f3 + transitions is 
independent of the value of LM parameter g' . Of course, 
however, the strength of the giant GT resonance state 
and each total strength of the /3_ and j3+ transitions de- 
pend on the value of the parameter. Let us estimate each 
strength separately. 



We divide the excitation-energy region into two 
parts: |g | < 2M* and |g | > 2M* . In the first region 
| go I < 2M*, we have a discrete state which corresponds 
to the giant GT resonance state in both NFA and NSA. 



In this region, we have 77p^ uli = 77Ao_ sca = 0, so that the 
RPA response function is described as 

7?RPA(go) = - — T-rlm 



(i) 



X 5 8tt 4 F T {q ) + ie 



with 



F T (qo) =Qo + X5 (- « P h + 9o77paiii(9o)) ■ 

In the case of NSA, 77 p ^ n (go) is replaced with 

77no- SC a(go)j but they are the same, as mentioned before. 
The excitation energy ujq of the giant GT resonance state 
is given by the solution of the equation: 

F T (u ) =0. 

When we expand 7 7 r(go) a t as, 

7 7 T(go) = 74(w )(go - w ) + O((g - ^o) 2 ), 

we can describe the response function in the form: 

1 V LOq 



7?RPA(go) 



X5 8^ 3 F^o 



■ S(q - w )- 



From this equation, we obtain the excitation strength of 
the GT state: 



S, 



1 V 



GT 



X5 8tt3 D' t (lj q ) ' 



(43) 
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d' t (u ) = X5 ( ^ + ^-iifiM) , and Sgt are a pp roximatel y § iven b y 

V 

using the dimesic function D T (qa) defined by uj q w xs^ph , Sgt ~ - 3 a p h- (44) 



(2tt) 3 



^t(?o) = <7o£>t(<7o)- 



Thus, the value of Sqt is almost independent of the value 
If we neglect n^ uU , which is small at |go| *C 2Af*, loq of the LM parameter. 

I 



In the region |go| > 2M* , the RPA response function in NFA is described as 



i?RP A (go) = ^ gsL&O f ffp(ftl) = _ rp + 

87r n _L^.n- ^„^ 2 _l /^.n-m r„_^ ®> 



(l+X5iIr(g )) 2 +(x 5 4au 1 i('?0) 

This equation provides us with the total GT strength of the /?_ transitions: 

/" V /' 2Sfp 

S NFA = / ^O-RrPAWo) = £-7 / d <70 — 

«» 2M* 87T J 2M , (1 



2BFp , /Pauli(fc p ,g0) 
^2M* ' "" 87T 4 7 2M , (1 + X5n r {qo)) + (X5^Pauli(fcp,go)) 

and that of the /3+ transitions: 



,( + ) /" , p , x V r 2 ^" , /Pa uli (fcn,g ) 

'NFA = / d 10 -KRPAl- ?0j - 5-4 / d 90-— — 2 " J 

•/ 2M* 87r J2Af* (1 + X5n r (-qoj) + (X5^Pauli(fen,go)) 



(45) 



where we have used the notations 
E Fi = 



M 



*2 



( i = p, n ). 



The bounds of integral are determined by Eq.()24p. 

The GT strengths of the region \q \ > 2M* in NSA 
are also calculated in the same way. They are obtained 
in terms of those in NFA as, 



c(-) _ „ ok 
J NSA ~ NFA ' 



?(+) 



(+) 

NSA 



(-) 
NFA- 



(47) 



Thus, each total strength of the /?_ and /3+ transitions 
in NSA is different from the one in NFA. We note that 
the above relationships are obtained by performing the 
integration of the NSA response functions over negative 
energy, as we did before. 

In Table |U we list the excitation energy and strength 
of the giant GT resonance state as a function of g' , and 
compare the strength of each energy region with the total 
RPA strength in NFA: 



total 



s, 



GT 



d(-) _ q(+) 
J NFA °NFA- 



In addition to (/tt/tott) 2 = 392MeV-fm 3 , we have em- 
ployed the values of parameters, as an example, M*/M — 
0.6, k n — 1.4fm _1 and k p = 1.2 fin . These give the 
sum rule value in the unit of the nuclear volume as 



N 



V 



— = ^{ k l- k l) = 0-06863 fnT 3 . (48) 

It is seen in the last column in Table [I] that this value is 
reproduced in NFA. Table"!] also shows that the strength 



TABLE I: Dependence of the excitation energy u>o and 
strength Sgt of the GT state on Landau-Migdal parameter 
g'. For details, see the text. 



9' 


Wo 


S G t/V 




s£L/v 


Stotal/y 




MeV 


fm" 3 


fm" 3 


fm" 3 


fm" 3 


0.0 


0.0 


0.06072 


-0.00733 


-0.01524 


0.06863 


0.5 


11.944 


0.06116 


-0.00750 


-0.01497 


0.06863 


1.0 


23.975 


0.06161 


-0.00767 


-0.01470 


0.06863 


2.0 


48.300 


0.06250 


-0.00804 


-0.01417 


0.06863 


3.0 


72.977 


0.06341 


-0.00844 


-0.01366 


0.06863 


4.0 


98.012 


0.06432 


-0.00887 


-0.01318 


0.06863 



of the giant GT resonance state does not depend strongly 
on the value of g 1 , as expected from the previous discus- 
sions. Comparing with sum rule value, Sgt is quenched 
by 12% ~ 10% for g' = ~ 1. From the experimental 
value for the excitation energy of the gian t GT resonance 
state, g' is estimated to be about 0.6[llj. In the present 
case, the approximate values given by Eg. 1(4*41 are 



UJQ 



23.801 g MeV, 



V 



0.06072 fm" 



9 



B. The Energy- Weighted GT Sum Rule 

In no charge-exchange excitations, there is a famous 
theorem on the energy- weighted sum value of the excita- 
tion strengths in RPA by Thoulessjij. According to the 
RPA theorem, the sum value is equal to the expectation 
value of the double commutator of the nuclear Hamilto- 
nian with the excitation operator. Here the expectation 
value is calculated by the ground state in the Hartree- 
Fock approximation. If the double commutator is a c- 
number, the sum value is given by the model-independent 
sum rule. Even if it is not a c-number, the theorem 
has been frequently used for a justification of the nu- 
merical results. In this subsection, we will show that 
the same theorem holds for charge-exchange excitations 
with respect to the sum value of the /3_ and /?+ transition 
strengths. 

In the case of NFA, the energy- weighted GT sum value 
is given by 



IV ( f°° 

- lim / dgogoIm#RPA(<3fo) 
it {Ziry r,^o \J_ V 

o 

dgogoImI?RPA(-<7o) 



7 (i) 

7 RPA 



where the first term of the r.h.s. stands for the energy- 
weighted sum of the /3_ transition strengths, while the 
second term that of the /3+ ones. In the same way as in 
the previous subsection, we expand /7rpa in terms of the 
coupling constant xs> 



(49) 



where In is defined as 



dq q Im.(n(qa)) 



n+l 



dq Q q Q lm (TT(-Q'o)) 



n+l 



(50) 



From the same arguments for Eq. (|40|) the above equation 
can be expressed as 



n+l 



/« = -Im / dq q (n(q )) n+1 + (iT(-g )) 
Jc 1 



At |3o| — > oo, the integration on the contour C gives 



dq q (II(±q )) 



n+l 



= r T2(2tt) 

i?™" 1 V V j jo 

Consequently, we obtain 

JW=0, (n>2). 



* N Z \ I d6e- 1 ^ 9 . (51) 



Thus the energy-weighted sum value is given by the two 
terms, 



1 V 



? (i) _ 1 V ( T (i) T (i)\ 



(2tt) 



where Iq 1 ^ is calculated using Ea . (|34|l . 



r(i) - 



while J-^ is obtained directly from Ea. (|51() . 

2 



-8tt d 3 p 



r-~ p p v (4(p) + a(n] 
p ~ u p 



r(i) 



7T 2(2tt) : 



■ N — Z 
V 



(52) 



(53) 



We note that Eq.JSTJ should not be used for Zq which 
requires the l/^o term in Ea. (|41l) . Finally, we obtain 
the energy-weighted sum of the /3_ and /?+ transition 
strengths in RPA of NFA, 



(1) 

RPA 



V 



d 3 p- 



)(n) 
P 



{N-Zf 
V 



(54) 



If the antinucleon degrees of freedom were neglected, 



we would have Iq — and I\ L> = — 7rap h , which give 



(i) 



*^RPA - 7^3^ 5a ph- 



r 



(2tt) 3 



(55) 



This is equal to uiqSgt from Eas. (|44|l . In fact, even if we 
take into account the coupling between the particle-hole 
states and the nucleon-antinucleon states, the energy- 
weighted strength of the giant GT resonance state is ap- 
proximately given by the above equation, 



V 



GT 



(27T) 



l - -vi 



(n - zy 

V 



(56) 



since the coupling is weak. This implies that the energy- 
weighted sum of the strengths in the giant GT reso- 
nance region is quenched by a factor (1 — 2w|/3) 2 , com- 
pared with the non-relativistic one. Its value is about 
0.77 for v F = 0.43 obtained from M* = 0.6M and 
k F = 1.36&0T 1 . 

It is worthwhile noting that we can do formally the 
above calculations including the divergent terms with 
GfGf- In this case, dp is replaced with 0p^ — 1 in 
Ea. (|52|) . so that we have 



S 



(i) _ 

RPA — 



V La P 



2 P 2 



v - (2-6^ 
E„ V p 



+4<7; 



(N — Z) 2 
V 



(57) 
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We can see that the negative contribution in Eq.JMJ) 
stems from neglect of the divergent term. A part of con- 
tribution from the Pauli blocking terms is included in the 
second term of Eq. l|54l) also, as seen from Eas. l|55|l and 

In NSA, we can perform the same calculation, but by 
defining instead of l£ } in Eq.® for NFA, 



negative spinors are given by 



J (1) - 



«<t(p) = 



E p + M* 



n+l 



dq Q q Q lm (H (q )) 
-n 

dq a q lm. (n(-q a )) 



2E 71 



E p + M* 



2E n 



1/2 



1/2 



1 

g P 

E p + M* 

v P 

E p + M* 
1 



ra+l 



with II(qo) given in Ea. i|35|) . We can show that the 
energy- weighted sum in NSA is the same as Ea. <|54|) for 
NFA. 

Now we will show that Eq. l|54|l is equal to the expec- 
tation value of the double commutator of the Hamilto- 
nian with the GT operator, when the expectation value 
is calculated by the ground state in the mean field. The 
nuclear field in the present model is written as 

V'(x) = J ( 2 ^3 /2 MP) exp(ip ■ x) a a (p) 

+v a (p) exp(-ip ■ x) bl(p)) , (58) 

where the first term in the parentheses describes the 
particles and the second the antiparticles. The suffix 
a denotes the spin and isospin quantum numbers as 
u a (p) = u a (p) |t) , ( a = a, t ), and the positive and 



J 



with the Pauli spinor £. The creation and annihilation 
operators, a a (p) and &J,(p), are defined as usual. Then, 
the mean field Hamiltonian, the interaction and the GT 
field operator are described, respectively, as 



H = d 3 x^(x)(~i-y ■ V + M*)V(x), 



V = 
F± = 



95 
2 



d 3 a;V(x)rf^(x)V(x)r Ml V(x) 



d 3 xip(x.)"/ 5 "f y T±ijj(x). 



We are assuming throughout the present paper that 
Dirac particles are bound in the Lorentz scalar and vec- 
tor potentials, but in the above mean field Hamiltonian, 
we have deleted the Lorentz vector potential which does 
not appear explicitly in the present discussions, while the 
Lorentz scalar potential is included in the nucleon effec- 
tive mass M*. 



The double commutator of the mean field Hamiltonian with the GT operator is easily calculated. By using the 
relationship for arbitrary operators A and B 1 



[V f (x)AV(x), V f (y)^(y)] =V f (x)[A, B]V(x)*(x-y), 



we obtain 



[F+, [ Hq , F- } ] = At 
ground sta 
\F + ,[Ho,F-]] |) 



i J d 3 x^(x)j ' V - j y -^\ V>(x) 



Its expectation value of the ground state in the mean field becomes 



(2-Kf 



J d 3 p [(e® + 9^) Tr CT (( 7 • p - lyPy ) A + (p)) 
+ 2Tr CT (( 7 -p- 7l/ ^)yl_(-p))] 



V 



2 2 

3 P -Py 



d J P 



E n 



2 - eg' 



(59) 



The double commutator of V with the GT operator is calculated in the same way. Using the abbreviations: 

r q = 70757% , XZ=[fg,fl], (q = ±,0), 

it is described as 



[F+, [V,F-]] 



95 

2 



d 3 x (^(X„ q )H tfXty + tff^ tf[fl,X» q ]il> 
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In the mean field approximation, the expectation value is calculated neglecting the exchange terms of the matrix 
elements. Keeping the only direct terms, we have 

< | [F+,[V,F-]] I > = -f Jd 3 x (( \^(X^ + )^^X^j\) + (\^X^j^(X^ + )^\ ) 

+ ( | tffai/, V f [ /+ , Xg ] ^ | ) + < | [ f + , Xg ] ^ | )) . 

The straightforward calculation of the above matrix elements yields 

(\[F + ,\V,F.}]\) = g b fcPxd^Hxm^Xl^l) 

2 = 4 95 (^!. (60, 



Thus, the sum of Eas. (|59|l and (|60|l is just equal to the 
energy-weighted sum of the strengths for the /?_ and /?+ 
transitions in RPA given in Eg. 15711 . 

sV A = {\[F + ,[Ho + V,F_}]\), 

including the divergent term. In Eq. l|54|l for NFA, the 
divergent term has been simply neglected. 

In NSA, the nuclear field is given by replacing the cre- 
ation operator of the antiparticles with the annihilation 
one in Ea. (|58|l for NFA. Because of this change, on the 
one hand, we obtain, instead of Eq.JSSJ), 

(\[F.,[H ,F + }]\) 

which does not contain the divergent term. On the other 
hand, the expectation value of the double commutator 
as for V is the same as Ea. (|60|) . Thus, in NFA also, the 
energy- weighted sum of the GT strengths in RPA is equal 
to the expectation value of the double commutator, when 
the expectation value is calculated with the ground state 
in the mean field. 

Formally we have proved that the RPA theorem holds 
in charge-exchange excitations also, but with respect to 
the sum of the strengths for the /3_ and f3+ transitions. 
In the present relativistic model, however, we have also 
shown that the energy-weighted sum value itself is diver- 
gent. If we neglect simply the divergent terms as in NFA 
and NSA, the sum value becomes negative, owing to the 
strengths of the Pauli blocking terms or the antinucleon- 
hole excitations. Generally speaking, all previous calcu- 
lations in NFA and NSA have the same problem. In order 
to solve this problem, we need definitely the renormaliza- 
tion of the divergence. 

V. DISCUSSIONS AND CONCLUSIONS 

The relativistic model has been extensively used as 
a phenomenological model of nuclei for the past 30 



years 0|. In particular, it explains very well the 
ground-state properties of nuclei with the mean field ap- 
proximation, where antinucleon degrees of freedom are 
neglected [la, Il3- In RPA based on the mean field ap- 
proximation, however, it is known that we can not de- 
scribe the excited states within the nucleon space only, 
and should include at least Pauli blocking terms from 
nucleon-antinucleon excitations in the correlation func- 
tions. We have called this approximation NFA. The Pauli 
blocking terms are required for RPA to keep the continu- 
ity equation 0, 0, and to reproduce the correct Landau- 
Migdal(LM) parameters, etc.Q. The reason why we need 
the Pauli blocking terms may be partially because of the 
fact that in relativistic models, the complete set needs 
antinucleon degrees of freedom. 

In this paper, we have shown that the Pauli blocking 
terms are also necessary for RPA to satisfy the Gamow- 
Teller(GT) sum rule with respect to the difference be- 
tween the strengths of the /3_ and /3+ transitions which 
is called Ikeda-Fujii-Fujita(IFF) sum rule [2 • This fact 
has been shown in an analytic way for nuclear matter. If 
the configuration space is limited to the nucleon one, the 
sum rule value is exhausted only by about 88%. When 
adding the Pauli blocking terms to the RPA correlation 
function, the sum rule value is reproduced. The coupling 
of the particle-hole states with the nucleon-antinucleon 
states is weak for reasonable values of the LM parameter 
g' . As a result, the GT strength which is distributed over 
the giant GT resonance region remains to be quenched 
by about 10 to 12%. 

In the previous paper [3j, the GT strength in the nu- 
cleon sector was estimated for finite nuclei in the mean 
field approximation. The value of the total strength was 
about 94% of the sum rule value. The reduction of the 
quenching is due to a larger value of the effective mass in 
finite nuclei than in nuclear matter. Since the coupling 
between the particle-hole and particle-antiparticle states 
is weak, the about 6% quenching of the GT strength is 
also expected in RPA of NFA for finite nuclei. 

The IFF sum rule in the no-sea approximation(NS A) 
has been also investigated. In NSA, the Dirac sea is 
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assumed to be empty, and the antinucleon states are 
treated as particle states with negative energy. In this 
way, one can avoid the divergence problem without vio- 
lating the continuity equation as in NFA. We have shown 
that NSA also satisfies IFF sum rule, and predicts the 
same quenching of the strength for the giant GT res- 
onance state. It has been shown, however, that each 
strength of the total /3_ and (3 + transitions is different 
from those in NFA. 

The energy- weighted sum of the GT strengths in RPA 
is also studied in NFA and NSA. It has been shown 
that the sum of the energy-weighted strengths for the 
j3- and /3+ transitions is equal to the expectation value 
of the double commutator of the Hamiltonian with the 
GT operator, when the expectation value is calculated 
with the ground state in the mean field and the divergent 
terms are deleted. Thus the well-known RPA theorem by 
Thouless^ holds for charge-exchange excitations also. 

We should note finally that renormalization of the di- 
vergence should be investigated in the future study of 
relativistic models. So far most of the nuclear observ- 
ables are well reproduced phenomenologically without 
the renormalization. NFA and NSA, which satisfy vari- 
ous conservation laws, are such examples. As shown in 
the present paper, however, NFA and NSA provide us 
with unphysical results also, as a price of neglecting the 
divergence. In the case of the GT excitations, the non- 
energy weighted strengths themselves in NFA and the 
energy-weighted strengths in NSA are negative outside 
of the giant resonance region. All previous calculations 



using NFA and NSA may have the same problems. 

It may depends on nuclear observables whether or not 
effects of the renormalization are important. In the previ- 
ous studies^, on the one hand, it was shown that the LM 
parameter F\ depends on antinucleon degrees of freedom, 
only through the Pauli blocking terms. These facts were 
shown by renormalized calculations in the a — lu model. 
As a result, some physical quantities which are domi- 
nated by Fi can be described using approximation with 
the Pauli blocking terms. Indeed, for example, as well 
as the center of mass motion, the space part of the nu- 
clear current, which is responsible for the orbital part of 
the nuclear magnetic moments was shown to be described 
well in NFA and NSA0,|22|. On the other hand, effects of 
the antinucleon degrees of freedom on the LM parameter 
Fq are not represented by the Pauli blocking terms only, 
but the contribution from other nucleon-antinucleon ex- 
citations is more important, as shown in ref . pi Il5j . An- 
other example is the Coulomb sum rule, where the renor- 
malization provides us with a strong quenching of the 
sum rule value at high momentum transfer|l5|. 

Thus, there may be some cases where the renormal- 
ization is not essential for description of the observables, 
but generally speaking, we should investigate in the fu- 
ture how effects of the renormalization change previous 
results in relativistic models. This may be also true for 
the GT strengths discussed here, although the coupling 
between the particle-hole terms and the Pauli blocking 
terms is weak in RPA. 
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